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Phase diagram and heavy-ion coUisions: Overview*^ 
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The physics of the QCD phase diagram is discussed in view from heavy-ion collisions, 
compact astrophysical phenomena, lattice QCD and chiral effective models. We find that 
(T, /x_b) region probed in heavy-ion collisions and the black hole formation processes covers 
most of the critical point locations predicted in recent lattice QCD Monte-Carlo simulations 
and chiral effective models of QCD. We also discuss the consistency of the statistical model 
results and dynamical model description. 

§1. Introduction 

While quarks and gluons are the fundamental degrees of freedom in the Quantum 
Chromodynamics (QCD), atomic nuclei are made of nucleons. When we heat-up 
or compress nuclei, the QCD phase transition takes place and quarks and gluons 
start to move over a large volume. Many of the current and past heavy-ion collision 
experiments have been carried out to discover the deconfined and thermalized matter 
made of quarks and gluons; the quark-gluon plasma (QGP). Recent experiments at 
the relativistic heavy-ion collider (RHIC) and the large hadron collider (LHC) have 
shown that strongly interacting matter with high energy density and characteristic 
transport properties is created,'D' and this matter is now commonly understood as 
a strongly coupled QCP-^^^ One of the important directions of current and future 
heavy-ion physics is to elucidate the properties of QGP and dynamics of heavy-ion 
collisions including its preequilibrium stage. Another intriguing direction is the study 
of the QCD phase diagram. 

QCD phase diagram has rich structure. At low temperature (T) and low baryo- 
chemical potential (/is), the hadron phase is realized where the chiral symmetry is 
spontaneously broken and color is confined. In the very large region, perturbative 
QCD predicts that the ground state is the color superconductor (CSC), where the di- 
quark pair condensates.® Recent large Nc arguments suggest the existence of another 
form of matter, confined high-density matter referred to as the quarkyonic matter.SJ 
These forms of matter can be probed in heavy-ion collisions and/or compact stars, 
as shown in Fig. [TJ 

The QCD phase transition has two facets, deconfinement and chiral transitions. 
The deconfinement transition involves a large degrees of freedom change. In the 
simplest bag model picture, pressures in the hadron gas and QGP are given as 
Ph = 3PsB and -Pqgp = 37PsB — B, respectively, where PsB = vr^T^/DO and B is the 
bag constant. The degrees of freedom suddenly increases from 3 (pions) to 37 (quarks 
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and gluons) at T = = (455/177r2)V4 ^ 0.725^/'^, which gives ~ 170 MeV 
for B^^^ = 240 MeV. The spontaneous breaking of the chiral symmetry and its 
restoration are understood as a consequence of the zero point energy and thermal 
contribution of quarks. We consider here a simple model for Nj = 2. The fermion 
contribution to the free energy density is given as, 

p. A'PA^ KPn^ NPT^ „, ,, , , 



dfV 



where M is the fermion mass, A is the cut-off, df = ANcNf, and Pj^ = 7/8 x 
TT^T^/gO + ^2^2/24 + ^^/487r2. The first term in Eq. ([FT]) represents the zero point 
energy, for each momentum k), and the second integral shows thermal con- 

tributions. Provided that the constituent quark mass is proportional to the chiral 
condensate a and the free energy density is given as f2/V = Qp/y + baCj'^ with a 
constant coefficient 60-1 we find that the second order chiral transition can take place 
at + /il/Svr^ = T"^ {fig = 3/i), which is shown by the dotted line in Fig. [2j These 
simple estimates roughly coincide with the critical temperature Tc obtained in the 
lattice QCD Monte-Carlo (MC) simulations and the chemical freeze-out boundary, 
as discussed later. The actual QCD phase transition has both of these natures and 
is known to be cross over for physical masses of two light {u, d) and one strange (s) 
quarks {Nf = 2 + 1^ as shown in Fig. [TJ 

In this proceedings, we discuss the QCD phase transition at high T and small 
baryon densities in Sec.[2l the critical point in SecO and the phase diagram structure 
of dense matter in Sec. [H We give a short summary in Sec. [5l 
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Fig. 2. Predictions of the critical temperature and critical point location in comparison with the 
chemical freeze-out points, measured photon temperature, and the swept region during the black 
hole formation. 



§2. Phase transition at small densities 

In high-energy heavy-ion cohisions, we can probe the transition from QGP to 
hadron phase in the high T and low fig region. At very high energy, gluons in 
the small x region governs the dynamics, and nuclei are regarded as the color-glass 
condensed (CGC) state.--' Thermalization takes place after a short time from the first 
contact,'^ and QGP is formed. Collective flows develop during the hydrodynamical 
expansion stagej^® and abundant hadrons are formed at the critical temperature 
T ^ Tc- Success of hydrodynamic^''^ and statistical modeld^ implies that local 
thermal equilibrium is achieved in the early stage at RHIC, and it enables us to 
discuss the phase transition in equilibrium. 

The cross over nature of the transition at small baryon density allows us to 
define Tc in several ways, such as the chiral susceptibility {xa)-, strange quark num- 
ber susceptibility (xs), the Polyakov loop (P), and the 0{N) scaling functions 
{0{N)). The transition temperature in the lattice MC simulation seems to be con- 
verging in the range Tc = (145-185) MeV. Around half of the uncertainty comes 
from the type of fermions. Simulations with staggered fermions give, for exam- 
ple, Tc{X(t) = 151(3)(3) MeV with stout-link improved staggered fermion^^^i and 
Tc{0{N)) = 157 ± 6 MeV with asqtad and HISQ/tree actions.l^ It is pointed out 
that the contamination of "heavy" pion at finite lattice spacing in the previous re- 
sult^^ would be the main source of discrepancy.'^ Other fermions (domain wall 
and Wilson fermions) result in Tc = (160-184) MeV.'^ I would like to mention 
that Tc extracted from the Polyakov loop seems to be systematically higher than 
the chiral transition temperature, Tc(xct)- For example Tc{P) = 176(3) (4) MeV is 
25 MeV higher than Tc{xct) with stout action.^'' This ordering supports the idea of 
deviated chiral-deconfinement transition boundary at finite fi^^ 
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It is interesting to find that the predicted Tc is close to the chemical freeze-out 
temperature in the statistical model fit, Tcfo ~ 160 MeVP^ The statistical model 
is based on the assumption that hadrons are produced in thermal and chemical 
equilibrium, and it describes the relative yields of normal hadrons very well only 
with two (T and /is) or three parameters. Approximate agreement of Tc and Tcfo 
is consistent with the elliptic flow calculation.'^ In a hydro+cascade hybrid model, 
the rapidity dependence of the elliptic flow {V2) is well explained when the switching 
temperature is chosen to be just below Tc^^^ Elastic scattering dominates in the 
hadronic cascade stage, and chemical composition would be kept approximately. 

It should be noted that we cannot determine Tc directly from Tcfo; ^cfo is 
the temperature in the hadron phase, and it should be smaller than Tc- In the 
symposium Gupta proposed a method to evaluate Tc based on data.l^ By comparing 
the lattice results with the RHIC-BES fluctuation data (variance, skewness and 
kurtosis) systematically, the critical temperature is evaluated as Tc = 1751^ MeV. 
Another important observation is the direct measurement of the QGP temperature. 
Hadrons interact strongly also in the hadronic stage, then we need direct photon 
and/or lepton spectra. PHENIX results of the direct photon spectrum show that 
significant excess is observed at pT < 3 GeV / c compared to scaled p + p results and 
the exponential fit in central collisions gives an apparent inverse slope parameter of 
Tphoton = 221 MeV.'^ The initial temperature is estimated to be Tinit = (300 — 
600) MeV depending on the thermalization time in hydrodynamical simulations. 
These temperatures are clearly higher than Tc estimated in lattice QCD and Tcfo 
evaluated in the statistical models, Tinit > Tphoton > Tc (or Tcfo)- This ordering 
provides another evidence of the formation of the QGP prior to hadronic freeze-out. 

§3. Critical Point Search 

While experiments at RHIC and LHC have been extremely successful in QGP 
hunting and revealing the properties of QGP, the transition at low baryon densities 
is cross over and we have not seen evidence of the real (first or second order) phase 
transition. The QCD phase diagram is expected to have a critical point (CP), which 
connects the cross over transition at low densities and the first order phase boundary 
at high densities as shown in Fig. [Tp^ Once the location of CP is identified, the 
global structure of the phase diagram is known. Thus the CP is regarded as the 
cornerstone of the QCD phase diagram. Experimental and observational inputs are 
essential for the CP hunting, since theoretical predictions of the CP location scatters 
in a wide region of the phase diagram as shown in Fig. [2pSli 

The lattice MC simulation is the ab initio framework to solve non-perturbative 
QCD, but it is difficult at present to obtain precise results at high densities because 
of the sign problem. The fermion determinant has the property T>*(/i) = D[—^) 
and can take a negative value at finite /i. This determinant is the importance MC 
sampling weight, and the cancellation of the weight leads to a large error. Exten- 
sive works have been done or are ongoing to overcome or to avoid the sign prob- 
l^^mmmmmmmmm while we have uncertainties from the sign prob- 
lem, we find a trend that lattice MC favors high CP temperature Tcp comparable 
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to the zero density critical temperature Tc{fi = 0). The multi-parameter reweighting 
methocPiJ predicts {Tcp,hcp) = (162 MeV,360 MeV)(LR04), where Tdn = 0) = 
164 MeV. Compared to their 2002 results, (Tcp, /icp) = (160 MeV, 725 MeV)(LR02), 
smaller quark mass and larger spatial lattice size result in smaller /xcp- The canon- 
ical ensemble method is applicable to a small size lattice, and a non-monotonic 
behavior oi fi/T would signal the first order transition. The CP location is esti- 
mated as (rcp,/icp) - (0.87rc,6r) (LC08^ and (rcp,;Ucp) = (0.927 T^, 2.60 Tc) 
(LC11),I23 while these results are obtained with heavier u, d quark masses than phys- 
ical values. The Taylor expansion m. ji/T combined with the reweighting predicts 
/iCP — 420 MeV (LT04).'22J On the other hand, the critical mass is found to decrease 
at small non-zero /i for Nf = 3, suggesting there is no chiral critical point.'^S 

Chiral effective models of QCD have been utilized to predict the phase dia- 
gram structure, including the existence of CP,'^ and they generally predict the 
CP location in lower T and larger /x region compared with lattice MC results. For 
example, NJL,ES Polyakov loop extended NJL (PNJL),E3 PNJL with 8 quark in- 
teraction (PNJLgjSSJ and Polyakov loop extended quark meson (PQM;P models 
predict (rcp,AiCp) = (22, 1095), (74, 1062), (129, 708) and (91,1005) (in MeV), re- 
spectively.'^ Here we have taken the vector-scalar coupling ratio of 0.2. 

In Fig. [51 we compare the predicted CP locations with the chemical freeze-out 
points at various incident energies and the (T, //g) region expected to be probed in 
heavy-ion collisions. We also show the region which would be probed during the 
black hole formation processes.® We find that most of the predicted CP locations 
are in the accessible region in heavy-ion collisions or black hole formation. LR02 and 
LC08 predictions are in the inaccessible region, but updated or more recent results 
(LR04 and LCll) are reachable. Next, a simple estimate of the phase boundary 
r2 + /il/Svr^ = r2 is found to roughly coincide with the predicted CP locations and 
chemical freeze-out line. We also find that the CP may be close to the chemical 
freeze-out line in -^/sj^ = (5 — 130) GeV range,!^ i.e. between the top AGS energy 
{Ei^c = 10.6 A GeV, y/s^ = 4.85 GeV) and the RHIC energy in the first run 
{y/sNN = 130 GeV), including the top SPS energy (-Einc = 158 A GeV , y/s^ = 
17.3 GeV). This expectation is consistent with collective flow studies; AGS energy 
heavy- ion collisions are well described in hadron-string cascade models ,'23 'ESI' 'EZJ and 
the elliptic flow at ^s^n > 130 GeV is successfully described in hydrodynamic^''^ 
rather than cascade.'^ 

The Beam Energy Scan (BES) program at RHICpS is promising, since the above 
incident energy range is fully covered. As discussed in this symposium,'32J the quark 
number scaling of the elliptic flow is found to hold at ^/s^ > 39 GeV, but not for cp 
meson at y/s^ = 11.5 GeV. This scaling appears from the collective dynamics in the 
partonic stage, then the above observation may suggest the onset of deconfinement 
in the range 11.5 GeV < ^/s^ < 39 GeV. RHIC-STAR collaboration also finds 
the reduction of the kurtosis compared with the hadron resonance gas model results 
at ^SjvAT ^ 20 GeV.I^ The distribution of the the order parameter X follows 
exp[— Seff (X)], where S^s is the effective action in which other variables than X are 
integrated out. In the vicinity of CP, Ses{X) becomes very flat around the value at 
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CP (Xc), and X will have a distribution exjp[—a{X — Xc)^]. This distribution lead to 
a negative kurtosis k, defined as the correlated part of the 4-th order correlation.'^SJ' 
This large fluctuation may be observed via baryon number fluctuations, since the 
chiral condensate couples with the quark number density at finite /i^ and the order 
parameter becomes a mixture of the chiral condensate and quark number density.!^ 
The sign of the kurtosis or the skewness for the proton number along the chemical 
freeze-out line is sensitive to the model details and still under debatePD' "We need to 
understand the behavior of these fluctuation observables both from theoretical and 
experimental sides. Theoretically, we need to go beyond the mean fleld treatment 
and to include fluctuation effects, and it is also important to understand how strong 
the baryon number couples with the order parameter. Experimentally, we need to 
obtain data with smaller error bars at smaller interval of incident energies. 

Non-monotonic behavior of several quantities is also observed at SPS; horn, step 
and dale^^ The horn structure is the sharply enhanced /tt^ ratio observed at 
-y/Sivjv — 8 GeV, and is proposed as the signal of QGP formation in statistical model 
based on the bag model The strangeness density is similar or even lower 

in QGP, provided that the initial A^A^ collisions make abundant strange hadrons. 
Then the /n'^ ratio will have a sharp peak at T ~ Tc, decrease during the mixed 
phase, and reach a constant value of QGP. The horn is also reproduced in a more 
recent statistical model analysis; Tcfo is found to saturate (step) but /xcFO continues 
to decrease at ^s^n — 8 GeV, then having u quark is disfavored at higher 
energies.'I^ Hadron-string cascade models cannot reproduce horn, but explain the 
step to some extent; a large part of particles are produced via string fragmentation 
at ^Sjviv ~ 5 GeV, then suppressed re scatterings during the formation time leads 
to the inverse slope saturation.fSSli'ESli xhe description of the horn is improved in the 
hadronic cascade-|-hydrodynamics hybrid model.'^Zl' The latter three suggest the need 
of faster thermalization than cascade, but QGP formation is not necessarily required. 
Further studies are necessary to pin down the onset of deconfinement. Beam energy 
and system size scan in NA61/Shine may be helpful for this purpose.^^ 

§4. Phase structure in dense matter 

Phase transition in cold dense matter is important in compact astrophysical 
phenomena. The baryon chemical potential in neutron star core is calculated to be 
/is ~ 1650 MeV {pb = 1.12 fm"^) for M^ax = 2.17 Mq star.H9 This chemical 
potential is much larger than the transition chemical potential to quark matter in 
many of the chiral effective models, fic = (1000 - 1110) MeY^ In super novae, 
it is suggested that the transition to quark matter after the collapse and bounce 
stage would lead to successful explosion and may be detected by the second shock 
peak in anti-neutrino.'^S' One of the recent interesting suggestions is the critical 
point sweep during the dynamical black hole formation processes .'32' The majority 
of massive stars (M > (20 — 25)Mq) may collapse to black holes without explosions 
(faint-supernova), where very hot (T ~ 90 MeV) and dense [pg ~ 5po) matter 
is formed.'^ Since the formed matter is isospin asymmetric, the isospin chemical 
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Fig. 3. First order phase transition surface in {T, jj,B , S/j.) space calculated with PQM is compared 
with the BH formation profile, at t = 0.5, 1.0, 1.344 sec after the bounce (double lines), t = 1.344 
sec corresponds to the time just before the black hole formation. 

potential 6fi = {^n — = (/^d — IJ-u)/'^ is finite and reduces Tqp- The black hole 

formation profile {T, jj,g, S/i) evaluated in i^-radiation hydrodynamic^^ is found to 
go through the critical point in asymmetric matter obtained in some of the chiral 
effective models as shown in Fig. [3l Thus if the critical point is in the low T and 
high fig region, it may be probed in the black hole formation processes. 

There are a variety of proposed phases in cold dense matter, such as nuclear 
superfluid, hyperon admixture, meson condensation, quarkyonic matter, inhomoge- 
neous chiral condensate, color superconductor, and baryon rich QGP. The variety of 
phases is related to the variety of interaction and condensate. The vector coupling,!^ 
determinant- type six fermi interactionjSSJ and eight fermi interactiorP^''^ would af- 
fect the phase structure at high densities. When the color superconductor joins the 
game, the pairing correlation plays a similar role to temperature in smoothing the 
fermi surface, and we may have another critical point at low tW^^^ 

In order to pin down the phase diagram structure in dense matter, the first 
principle approaches are of course desired. One of the promising directions is the 
strong coupling lattice QCD,I2S1' where we expand the lattice QCD effective poten- 
tial in 1/(7^. Recently, the phase diagram is obtained with finite coupling effects 
(1/5^, in the mean field approximatiorpS' and with fluctuation effects in the 

strong coupling limit .'^S Both of these results are obtained in the strong coupling 
region (/3g = 2Nc/g'^ < 4 and /3g = 0) and based on the unrooted staggered fermion 
corresponding to Nf = 4 in the continuum limit. These are still far from realistic, 
but combination of these techniques may be helpful to understand the whole struc- 
ture of the QCD phase diagram. Another promising direction may be to improve 
effective models by using the functional renormalization group approaches. ^SSJ 



In this proceedings, the physics of the QCD phase diagram is reviewed with 
emphasis on the relation to heavy- ion collisions. Lattice MC simulations are powerful 
for high T and low fi transition, and may be useful to discuss the critical point (CP). 



§5. Summary 
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The CP location in chiral effective models depends on the model details, but it 
roughly follows a simple estimate based on the quark zero point energy and thermal 
contribution. Both hydrodynamics and statistical models are successful at RHIC, 
and we can assume local thermalization. These theoretical estimates of Tc and CP 
are compared with statistical model results and the {T,^g) region probed in heavy- 
ion collisions and black hole formation processes. If we are not very unlucky, CP 
would be in the accessible region in either of two. 

As pointed out in the symposium, thermalization may not be achieved at SPS or 
lower energies. Experimentally, it is necessary to cross-check the results depending on 
statistical model arguments by using more dynamical observables. Theoretically, it 
is challenging to develop theoretical frameworks which can describe non-equilibrium 
transport properties and the change of relevant degrees of freedom across the phase 
boundary. I could not mention on the chiral magnetic effect and exotic hadron 
production in heavy-ion collisions in this proceedings. 

Acknowledgement 

The author would like to thank Prof. A. Nakamura for useful discussions. This 
work is supported in part by the Global COE Program "The Next Generation of 
Physics, Spun from Universality and Emergence", and the Yukawa International 
Program for Quark-hadron Sciences (YIPQS). 

References 

1) K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757 (2005), 184; J. Adams et 
al. [STAR Collaboration], Nucl. Phys. A 757 (2005), 102; K. Aamodt et al. [The ALICE 
Collaboration], Phys. Rev. Lett. 105 (2010) 252302. 

2) M. Gyulassy and L. McLerran, Nucl. Phys. A 750 (2005), 30. 

3) M. G. Alford, A. Schmitt, K. Rajagopal and T. Schafer, Rev. Mod. Phys. 80 (2008) 1455. 

4) L. McLerran and R. D. Pisarski, Nucl. Phys. A 796 (2007) 83; 

5) Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, Nature 443 (2006) 675. 

6) L. D. McLerran and R. Venugopalan, Phys. Rev. D 49 (1994) 2233; E. lancu, K. Itakura 
and L. McLerran, Nucl. Phys. A 708 (2002) 327. 

7) U. W. Heinz and P. F. Kolb, Nucl. Phys. A 702 (2002) 269; Z. Xu and C. Greiner, 
Phys. Rev. C 71 (2005) 064901; A. Dumitru and Y. Nara, Phys. Lett. B 621 (2005) 
89; R. J. Fries, B. MuUer and A. Schafer, Phys. Rev. C 79 (2009) 034904; T. Kunihiro, 
B. MuUer, A. Ohnishi, A. Schafer, T. T. Takahashi and A. Yamamoto, Phys. Rev. D 82 
(2010) 114015. 

8) P. Huovinen et al. Phys. Lett. B 503 (2001) 58; 

9) T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey, Y. Nara, Phys. Lett. B636 (2006) 299. 

10) A. Andronic, P. Braun-Munzinger and J. Stachel, Nucl. Phys. A 772 (2006) 167. 

11) Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643 (2006) 46. 

12) A. Bazavov et al. [HotQCD Collaboration], J. Phys. G 38 (2011) 124099. 

13) M. Cheng et al, Phys. Rev. D 77 (2008) 014511. 

14) K. Kanaya, PoS LATTICE 2010 (2010) 012. 

15) M. Cheng et al., Phys. Rev. D 81 (2010) 054510; V. G. Bornyakov et al., Phys. Rev. D 
82 (2010) 014504; S. Ejiri et al. [WHOT-QCD Collaboration], Phys. Rev. D 82 (2010) 
014508. 

16) K. Miura, T. Z. Nakano, A. Ohnishi and N. Kawamoto, arXiv:1106.1219 [hep-lat]. 

17) S. Gupta et al, Science 332 (2011) 1525; S. Gupta, in this proceedings. 

18) A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 104 (2010) 132301; Y. Yam- 
aguchi, in this proceedings. 



Phase diagram and heavy-ion collisions: Overview 



9 



19) M. Asakawa and K. Yazaki, Nucl. Phys. A 504 (1989) 668. 

20) M. A. Stephanov, Prog. Theor. Phys. Suppl. 153 (2004) 139. 

21) Z. Fodor, S. D. Katz, JHEP 0203 (2002) 014. JHEP 0404 (2004) 050. 

22) S. Ejiri et al, Prog. Theor. Phys. Suppl. 153 (2004) 118. 

23) P. de Forcrand and O. Phihpsen, Nucl. Phys. B 642 (2002) 290; M. D'Eha and M. - 
P. Lombardo, Phys. Rev. D 67 (2003) 014505; 

24) S. Ejiri, Phys. Rev. D 78 (2008) 074507. 

25) A. Li, A. Alexandru and K. -F. Liu, Phys. Rev. D 84 (2011) 071503. 

26) K. Nagata and A. Nakamura, Phys. Rev. D 82 (2010) 094027. 

27) P. de Forcrand and O. Philipsen, JHEP 0701 (2007) 077; JHEP 0811 (2008) 012. 

28) T. Z. Nakano, K. Miura and A. Ohnishi, Prog. Theor. Phys. 123 (2010) 825; Phys. Rev. 
D 83 (2011) 016014; K. Miura, T. Z. Nakano and A. Ohnishi, Prog. Theor. Phys. 122 
(2009) 1045; K. Miura, T. Z. Nakano, A. Ohnishi and N. Kawamoto, Phys. Rev. D 80 
(2009) 074034. 

29) P. de Forcrand and M. Fromm, Phys. Rev. Lett. 104 (2010) 112005; W. linger and P. de 
Forcrand, J. Phys. G 38 (2011) 124190. 

30) Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; Phys. Rev. 124 (1961) 246; 
T. Hatsuda and T. Kunihiro, Phys. Rept. 247 (1994) 221. 

31) K. Fukushima, Phys. Lett. B 591 (2004) 277; S. Roessner, C. Ratti and W. Weise, Phys. 
Rev. D 75 (2007) 034007. 

32) K. Kashiwa, H. Kouno, M. Matsuzaki and M. Yahiro, Phys. Lett. B 662 (2008) 26; 

33) B.-J. Schaefer, J. M. Pawlowski and J. Wambach, Phys. Rev. D 76 (2007) 074023; 
V. Skokov, B. Friman, E. Nakano, K. Redlich and B.-J. Schaefer, Phys. Rev. D 82 (2010) 
034029. 

34) A. Ohnishi, H. Ueda, T. Z. Nakano, M. Ruggieri and K. Sumiyoshi, Pliys. Lett. B 704 
(2011) 284. 

35) S. A. Bass et al., Prog. Part. Nucl. Phys. 41 (1998) 255; Y. Nara, N. Otuka, A. Ohnishi, 
K. Niita and S. Chiba, Phys. Rev. C 61 (2000) 024901; P. K. Sahu, W. Gassing, U. Mosel 
and A. Ohnishi, Nucl. Phys. A 672 (2000) 376; P. Danielewicz, R. Lacey and W. G. Lynch, 
Science 298 (2002) 1592; E. L. Bratkovskaya et al, Phys. Rev. G 69 (2004) 054907; M. Isse, 
A. Ohnishi, N. Otuka, P. K. Sahu and Y. Nara, Phys. Rev. G 72 (2005) 064908. 

36) P. K. Sahu, S. G. Phatak, A. Ohnishi, M. Isse and N. Otuka, Pramana 67 (2006) 257. 

37) M. Nahrgang, G. Herold, S. Schramm and M. Bleicher, arXiv:1103.0753 [hep-ph]. 

38) T. Hirano, M. Isse, Y. Nara, A. Ohnishi and K. Yoshino, Phys. Rev. G 72 (2005) 041901. 

39) G. Odyniec, Acta Phys. Polon. B 40 (2009) 1237; D. Gebra, in this proceedings. 

40) H. Fuju, Phys. Rev. D 67 (2003) 094018. 

41) M. Asakawa, S. Ejiri and M. Kitazawa, Phys. Rev. Lett. 103 (2009) 262301; B. Friman, 
F. Karsch, K. Redlich and V. Skokov, Eur. Phys. J. G 71 (2011) 1694; M. A. Stephanov, 
Phys. Rev. Lett. 107 (2011) 052301. 

42) S. V. Afanasiev et al. [The NA49 GoUaboration] , Phys. Rev. G 66 (2002) 054902; G. 
Stefanek, in this proceedings. 

43) M. Gazdzicki and M. I. Gorenstein, Acta Phys. Polon. B 30 (1999) 2705. 

44) T. Gzopowicz, in this proceedings. 

45) G. Ishizuka, A. Ohnishi, K. Tsubakihara, K. Sumiyoshi and S. Yamada, J. Phys. G 35 
(2008) 085201. 

46) T. Hatsuda, Mod. Phys. Lett. A 2 (1987) 805. I. Sagert et al. Phys. Rev. Lett. 102 (2009) 
081101. 

47) K. Sumiyoshi, S. Yamada, H. Suzuki and S. Ghiba, Phys. Rev. Lett. 97 (2006) 091101. 
K. Sumiyoshi et al. Astro phys. J. Lett. 690 (2009) L43; K. 'i. Nakazato et al, Astrophys. 
J., to appear [arXiv:1111 .2900 [astro-ph.HE]]. 

48) M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, Prog. Theor. Phys. 108 (2002) 929. 

49) M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 44, 1422 (1970); G. 't Hooft, Phys. 
Rev. D 14 (1976) 3432 [Erratum-ibid. D 18 (1978) 2199]. 

50) B. Hiller et al, Phys. Rev. D 81 (2010) 116005; B. Killer, in this proceedings. 

51) Z. Zhang, K. Fukushima and T. Kunihiro, Phys. Rev. D 79 (2009) 014004; T. Kunihiro, 
in this proceedings. 

52) T. K. Herbst, J. M. Pawlowski and B. -J. Schaefer, Phys. Lett. B 696 (2011) 58; K. Morita, 
V. Skokov, B. Friman and K. Redhch, Phys. Rev. D 84 (2011) 074020. 



